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Abstract 

> : 

! Magnetic fields (MF) in the Local Supercluster (LSC) of galaxies may have 

^ I profound consequences for the propagation of Ultra High Energy Cosmic Rays 

OO I (UHECR). Faraday rotations measurements provide some informations about MF 

^ ' in compact clusters. However, very few is known about less dense regions and 

O ■ about the global structure of MF in the LSC. In order to get a better knowledge of 

^ ■ these fields we are performing constrained magnetohydrodynamical simulations of 

the LSC magnetic field. We will present the results of our simulation and discuss 

P ! their implications for the angular distribution of expected UHECR deflections. 

+3 . 
c/3 . 

• 1. Introduction 

^> : 

I Several arguments suggest that Ultra High Energy Cosmic Rays (UHECR) 

I are extragalactic in origin (see e.g.[l]). Propagation of electrically charged CR 

■ - - ' primaries is affected by intergalactic MF. Resulting picture crucially depends on 

magnetic field mean intensity and power spectrum. MF have been observed in 
several clusters of galaxies [7]. An indisputable evidence of these fields is provided 
by diffuse cluster-wide synchrotron emission from several clusters compatible with 
a field strength larger than 0.1 fiG. In addition, indirect observational evidence 
of the magnetization of the Inter Cluster Medium (ICM) comes from Faraday 
Rotation Measurements (RM) of the polarized radio sources located within or 
behind the clusters. Analyzing these observations several authors (see e.g. [3]) 
found that the ICM in clusters are permeated with a high filling factor by magnetic 
fields at a level of 1 — 10 fiG with a correlation length of 10-^ 100 kpc extending up 
to 1 Mpc from the cluster center. Beyond this distance the IGM density becomes 
too low to allow an observable Faraday rotation since the effect is proportional 
to the free electron density. In this case only upper bounds are available which 
are B < 10~^ 10~^ G in unclustered regions and a significantly weaker upper 
limit B < 10~^ in the sheets or the filaments of galaxy clusters [2]. What is most 
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crucial for UHECR reaching the Earth is the MF in the Local Supercluster (LSC) 
since it comprise most of the GZK volume. The MF structure of the LSC is very 
poorly known. RM allowed to determine its intensity only in some overdense 
regions like the Virgo cluster where it was found to be about 10~^ G. A method 
to determine the power spectrum of IGMF from a systematic analysis of RM of 
a large number of distance sources may allow a better knowledge of the LSCMF 
in a not too distant future [5]. Meanwhile, however, the only way to improve 
our understanding of UHECR observational data is to estimate the LSCMF large 
scale structure by numerical simulations. 

2. MHD simulations in clusters 

Magneto Hydrodynamical (MHD) simulations of MF evolution in galaxy 
clusters have been already performed by one of us (KD) and collaborators [4]. 
This kind of codes combines the merely gravitational interaction of a dominant 
dark-matter (DM) component with the hydrodynamics of a magnetized gas to 
simulate the formation of magnetized galaxy clusters. Initial conditions are spec- 
ified at redshift Zin ~ 20 by a random choice of density perturbations compatible 
with standard ACDM cosmology and an initial seed field to be tuned to repro- 
duce RM of clusters at redshift close to zero. The MF is amplified by adiabatic 
compression during cluster collapse and, more efficiently, by the Kelvin-Helmholtz 
instability where merger events give rise to strong shear flows. It was found that 
the simulation reproduce quite well the general features of the observed RM with 
an initial fleld strength of 5 x 10~^ G. In the outer regions of the clusters, the av- 
erage fleld strength follows the gas density proflle quite well though with a steeper 
slope compared to estimates based on flux conservation and adiabatic compres- 
sion (B oc Pg^s), showing that shear flow have to play an important role there. As 
a consequence the required initial intensity of the MF is smaller than expectation 
based on a mere adiabatic ampliflcation. The MHD simulations are designed to 
study individual clusters within a cosmological environment, and therefore can 
not describe voids or fllaments. Nevertheless an order of magnitude estimate can 
be obtained on the basis of the following simple argument. In the voids, as well 
as in low density regions around the LG, we do expect the adiabatic compression 
to give a rather good description of the MF evolution. Using flux conservation, 
which implies B{z = 0) = i?(2;in)(l + 2^111)"^, we see that a fleld strength of 5 x 10~^ 
G at Zin = 20, as required for the MHD simulations to reproduce RM of clusters, 
corresponds to the present time intensity of B{z = 0) ~ 10~^^ G in the unclus- 
terized IGM. The rms deflection angle of a particle of charge Ze traveling over 
the distance d through an irregular MF with a coherence length Ic <^ d and local 
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intensity B is given by the expression 

It is clear that a MF with strength ~ 10^^^ G have a neghgible influence on 
UHECR propagation. Therefore we expect significant deflections only when 
UHECR trajectory crosses cluster or a strongly magnetized filament. In that 
case deflection of few degrees can be produced (see FigH)). A different situation 
was examined recently by Sigl at al. [9] . The authors of [9] performed an Eulerian 
hydro+N-body simulation of the entire LSC and added to it a MF which follows 
passively the gas. Such simulation, however, is not a constrained simulation (see 
below) which implies that it is not suited to reproduce the real spatial distri- 
bution of matter and of magnetic fields. In [9] the authors chosen the observer 
position to lie in a region with a rather strong magnetization B ~ 0.1/xG. This 
was required in order to reproduce the observed isotropy in the arrival directions 
of UHECR under the hypothesis that the source distribution trace that of the 
matter in the LSC. Such an hypothesis, however, has not been proved so far and 
other reasonable possibilities have been discussed in the literature. In our opin- 
ion the assumption of a widely extended MF with intensity around 10"'' G in the 
surrounding of the LG is unrealistic. Indeed the LG lie in a quite peripherical 
region of the LSC and it is neither within a rich cluster nor within a strongly radio 
emitting region. So far these are the only extended systems where observational 
evidence of MF with strength comparable or larger than 0.1 yuG have been found. 

3. Constrained simulation in the LSC 

Realistic maps of expected UHECR deflections by EGMF as a function 
of distance to plasuable sources are highly desirable. Ideally, such maps should 
reflect obsreved distribution of galaxies and clusters and therefore different kind of 
simulations is needed. A promising approach is that of constrained simulations, 
i.e. simulations with initial conditions constrained to reproduce the observed 
large-scale structures in the nearby universe. Simulations of this kind have been 
already performed for the coUisionless DM and the gaseous components of the 
LSC (see e.g. [6,8]) which succeed to recover the main observed structures starting 
from initial conditions at high redshift compatible with standard cosmology. 

We are now working to produce a fully MHD constrained simulation of 
the LSC. Our code incorporates the MSPH technique [4], which was developed 
to follow MF evolution in simulated galaxy clusters, into a constrained N-body 
simulation of the DM and the gas in the LSC. Hopefully, at the end we will have 
simulated 3D maps resembling real local structure of MF fields. Once we will be in 
possession of such maps we will trace UHECR trajectories through the magnetic 
web. We will not limit ourselves by UHECRs sources located inside the LSC. 
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Deflection angle vs direction on the sky | 
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Fig. 1. Deflection angle as a function of the direction. 



In particular, we are interested to find out if and when correlations of charged 
UHECR with BL Lacs, found in Ref. [10], are compartible with expected EGMF 
structure. In other words, one of our main goals is to produce maps of expected 
UHECR deflections for a wide range of source distances. 

We will present detailed results of our simulation in the course of ICRC2003. 
We provide here a sample of our simulation. Fig. 1 represents a map of deflections 
of protons with energy E = 4 x 10^^ eV derived from a Virgo like cluster, seen 
from 5 Mpc. The MP intensity at the cluster center corresponds to ~ 1 /^G which 
is representative of a clusterized region. Rays have been propagated about 10-15 
Mpc in each direction. Deflections larger than 5 degrees are cut at this value to 
increase dynamical range of the plot for small deflections. 
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